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The properties of surface Dirac Fermions on a 3D topological insulator in proximity to a magnetic 
insulator with spatially textured magnetization are considered. We present an exact analytical 
treatment of the spectrum, the bound states and the domain wall resistance for an extended generic 
domain wall with in-plane and out-of-plane magnetizations. In the latter case, we find oscillations 
in the domain wall resistance as a function of the wall width and for certain widths a complete 
absence of reflections for all incoming momenta. The surprising occurrence of oscillations and 
the reflectionless potentials can be related to a supersymmetry of the surface Dirac Hamiltonian 
combined with the domain wall proflle. 
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Three dimensional strong topological insulators (3D- 
STI) are a recently discovered class of materials, with a 
prominent example being Bi2Se3 [Tj. Contrary to ordi- 
nary insulators, they exhibit topologically protected sur- 
face states with characteristic spin-momentum coupling, 
a result of strong spin-orbit interactions [2]. To exploit 
the full potential of these materials and for possible ap- 
plications, a combination of the surface states with more 
conventional materials like ferromagnets or superconduc- 
tors in proximity structure are desired. E.g., induced 
superconductivity gives rise to Majorana Fermions 
and induced magnetization textures exhibit a quantized 
magneto-electric effect [71 [8] . Ref. |9| is a recent review 
covering topological states of matter. 

Proximity induced ferromagnetism, where the order 
parameter can in general be inhomogeneous and time- 
dependent gives rise to phenomena interesting for spin- 
tronics and magnetotransport [lOlfTS] . Due to the spin- 
momentum locking, electrical current flow leads to a sig- 
nificant contribution to the spin-torque acting on the 
magnetization dynamics pTkl3) . Magnetically doped 
3D-STI could be used as a condensed matter realiza- 
tion of axion-electrodynamics [Hj. The transport of 
Dirac-Fermions through DWs has been also studied in 
Graphene [TS]. 

In this Letter, we consider a transport configuration 
in which a single static domain wall (DW) is located be- 
tween two contacts as illustrated in Fig. [T] We calculate 
the ballistic conductance for in-plane (IP) and out-of- 
plane (OOP) wall configurations. As our main result, 
in the OOP case, we find characteristic oscillations as 
a function of the wall width and/or the strength of the 
induced exchange potential. Such a signature for the 
surface states in transport experiments could be visible, 
even if residual bulk carrier density transport is present 
[16]. Interestingly, we find that the DW constitutes a 
reflectionless potential for certain wall widths. 

The general setup we study is illustrated in Fig. [T] 
and described by the effective Hamiltonian for the two- 




FIG. 1. The experimental setup for studying ballistic trans- 
port through a domain wall. Two leads with external voltage 
bias V are attached to the structure of topological insula- 
tor coated with a thin ferromagnetic layer that contains an 
in-plane (IP) tail-to-tail DW or out-of-plane (OOP) DW. 

dimensional surface electrons 

H = ihv(T-{e^xV) - Mm{r)-(T . (1) 

Here, the first term is the dispersion of the surface Dirac 
states |17| and the second term is the proximity induced 
exchange coupling to the magnetization profile m{r) 
with constant magnitude M . Such a magnetization tex- 
ture will occur naturally in films of magnetic materials, 
which we assume to be placed on top the topological in- 
sulator. The structure can in some limits be manipulated 
by an external magnetic field. We stress that the induced 
magnetization affects only the surface states and leaves 
the bulk conductivity of the 3D-STI unaffected - thus 
opening a possible path for disentangling the surface and 
bulk contributions to the conductivity. 

The in-Plane Configuration The lower domain wall 
sketched in Fig. [l] has the explicit form m{r) — 
(cosi?(a;), sin79(a::), 0). We assume the angle d{x) has the 
analytical form cos'd{x) = tanh (x/w). This shape can 
be obtained within a mean-field model with ferromag- 
netic exchange constant J and anisotropy constant K, so 
that the length of the domain wall becomes w = \J J jK 
|18j . Our problem is effectively one- dimensional and ky 
is a good quantum number. We make the ansatz for the 
wave-function |^'(a;,2/)) = e^fi^^'^^^e*^'^* 1-0(3;)), by which 
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FIG. 2. (Color online) Energy spectrum of the system with 
in-plane (a) and out-of-plane DW (b) including the bound 
states for Mw = i.'dhv. The spectrum is symmetric with 
respect to positive and negative energies in both cases. The 
continuum of scattering states (yellow area) has an energy 
gap 2M in both cases. Continuum states that exist only on 
either side of the DW (red and green areas in (a)) are decaying 
into the potential step posed by the DW. The inset illustrates 
various momenta relevant for the scattering states. 



the mj,-component is eliminated by the gauge factor, and 
end up with 



IP 



ihvaydx + {hvky — Mta,n\\{x/w))iJx ■ (2) 



From now on, we choose units such that hv — 1 and 
restore them only in the final results. 

In the homogeneous parts, we observe that the mag- 
netization shifts the Dirac cone along the fcj^-direction, 
which is also evident in the spectrum of this system 
shown in Figure [2^. The red and green cones correspond 
to the dispersions far away from the DW, so that the wave 
vectors in transport direction for the left (L) and the right 
(R) side obey fc^/R = — {ky ± M)^ when we consider 
an eigenstate with energy E. The current operator is 
j = {jxijy) = v{~ay, a^)- We find for incoming plane- 
wave states that the current points along (fcL/R, ky ± M) 
corresponding to an angle 7^^^^ = arctan((A:y ± M)/ki^/j,) 
(see inset Fig. [2^). 

We can solve the full eigenproblem Hip\4>) = E\iJ;) 
analytically by the substitution z = ^(l — tanh;^), 
which gives two decoupled difi'erential equations for = 
(^(fi^ (fil) . In order to bring these equations into hy- 
pergeometric form, we make the ansatz ip^{z) = (1 — 



(z) and find that 



{z(l - z)d^^ + (7 - (a + /? + l)z)d^ -a/3}(j>^0 (3) 



l + (M-iko)w, (3 = -(Af + 
using ko = ^{^L + ^r)- This 



with the coefficients a = 
ik())w and 7 = 1 — ikj^w 
equation is known 19J and we directly write the solution 



2 2 ('^R 2 ) 



(1 



2Fi(a,/3,7,z) , 
(4) 



where 2F1 is the hypergeometric function. The solution 
for (f^ can be obtained directly by using the symmetry 
[TnTkyCry^Hip] = 0, so that ipi{z) = -iTnTkyf^iz) (we 
define symmetry operations as Tm ■ M — >■ —M and simi- 
lar for ky and x). The prefactor of this solution is chosen 
such that it directly reflects the symmetry properties of 
the system, i.e. rMTk^ayl'il^'^^lj = 

Eq. ([4| exhibits the different types of solutions illus- 
trated in Fig. [2^, depending on whether /cl/r are real 
or imaginary. In particular, the bound state spectrum 
follows for fct and fc^ being imaginary and we straight- 
forwardly find 2[AJ -I- 1 bound states at energy i?„^o — 
±!i^J(2A-n)n, where we defined A = ^ ^ ^ 
with the magnetic length Im- The dispersion along ky is 

En.ky = En,oJl — w^k^/{A — nY and the states merge 



with the continuum at k. 



in) 

y,max 



(A 



7(wA). Here, 



n = describes a zero energy flat band for \ky\ < M as 
also seen in Fig. [2^. Note that the square-root depen- 
dence of the energy on the quantum number n is typical 
for Dirac Fermions, as for example in the case of Landau 
levels for Dirac Fermions [20] . 

The scattering solutions are described by real fcL/n 
and are pairwise degenerate. We can use parity sym- 
metry in the x-y-plane, i.e. [i?ip, Tfej^T^cTz] = 0, in order 

to obtain the second orthogonal solution, |V'_fep^ k ) ~ 

~'i-Tky%'^z\tpi^l )• Here, 1^'^^)^ ) describes an incoming 
wave from the left that is partialfy reflected and transmit- 



ted to the right side, and likewise, \ip. 



(IP) 

-fcR,fe„ 



describes an 



incoming wave from the right side. From the asymptotic 
behavior of Eq. (|4|, we can extract the transmission and 
reflection amplitudes and find the transmission probabil- 
ity 



Tip{E,ky 



sinh(7r/cLw) sinh(7rfcpi w) 
sin2(7rA) + sinh^d (fcL + k^)w) 



(5) 



and the refiection probability Rip = 1 — Tip. In the 
regime of small externally applied voltage V , we calcu- 
late the linear conductance G using the Landauer for- 
mula G = Gq ^(^F, ^y), where Gg = ^ is the 
fundamental conductance quantum [21j . In the absence 
of the domain wall, the conductance for transport along 
the x-direction is Go = Gq and W is the transverse 
dimension of the ballistic contact. 

To quantify the change of the conductance due to the 
presence of the domain wall, we define the domain wall 
resistance as 



SG^- 



G 



DW 



Gf) 



Gn 



= SG 



M 



SG 



DW ■ 



(6) 



We split this into two contributions: (5Gdw depends on 
the specific domain wall profile while 6Gm — ^ > ^ 
is the fraction of totally reflecting channels (red shaded 
area in Fig. [2^) to the total number of transport chan- 
nels (yellow + red areas). Essentially, 5Gm encodes the 
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FIG. 3. Relative decrease in the conductance 5G due to the 
presence of the domain wall for in- and out-of-plane DWs. In 
the case of an in-plane DW, we see that the values are rather 
large even for wider walls and decrease monotonically towards 
the asymptotic value 5Gm as the wall width increases. To 
the contrary, for out-of-plane DWs it decreases significantly 
when the wall width becomes of the order of the magnetic 
length Im ~ hv/M. The approximate formula (10 1 is drawn 



as black thin lines and shows very good agreement with the 
exact integration. 



spectrum mismatch between both sides of the wall, as 
far as it concerns transport at the Fermi level. The re- 
maining contribution from the DW is then related to the 
reflection from the wall profile and is explicitly given 
by (5Gdw = J-Ief-m)'^^ RipiEy.ky). The total 
change in conductance 5G is shown in Fig. [3j and we 
find that 5Gm strongly dominates, except for short walls 
for which w < Im- As the Fermi level moves deeper into 
the metallic regime, we also see an overall decrease of the 
change in conductance. 

The limiting behaviors can be obtained analytically. 
For sharp walls, i.e. w <^ Im, we can approximately 
solve the integral and find (5Gdw ""^^^ |(5Gm. Since 
SGuw decreases when the wall width increases, we can 
conclude (5Gdw ^ j^Gm- In the limit of wide walls, i.e. 
w » Im, using the saddle point method, we find 



2,, 2 



{2TTEFy 



1 



(7) 



The ballistic domain wall resistance 5Guw decays with 
the inverse square of the domain wall width, but the con- 
tribution SGm — is also present in this limit and is 
the dominating one. Experimentally, one could subtract 
SGm (for example determined in the absence of the DW) 
in order to obtain the contribution from the domain wall 
profile. 

As the Fermi level comes close to the band edge Ep « 
M, we have 6G — >■ 1, i.e. the domain wall blocks all 
transport channels. However, when the Fermi level Ep 
approaches M and eventually reaches the crossing point 
of the left and right Dirac cones (crossing of red and 



green line in Fig. [2p,), the number of states available for 
transport drastically reduces and one reaches the point 
of minimal conductivity. At this point, transport is due 
to few evanescent modes [HI [53], which is beyond the 
present study. 

The out-of-Plane Configuration. The OOP wall has 
the magnetization profile m(r) = (0, sini9(a;), cos-!^(x)) 
which, in the language of the 2D-Dirac equation, de- 
scribes a mass domain wall connecting two quantum 
anomalous Hall states of opposite chirality. We per- 
form a spin rotation around the y-axis by ^ utilizing 
the unitary spin rotation matrix U = e^^T'^y to obtain 
the representation U HqopU^ — Hip{ky = 0) — hvkyO-z, 
which allows us to reuse the previous results. As in- 
ferred from the symmetry {iJip,^^} = 0, we see that 
fcyCTz only couples pairs of positive and negative en- 
ergy. Thus, we only have to diagonalize 2x2 sub-blocks 
and straightforwardly obtain the full energy dispersion 



±y + h?v^{k'^ + k^) with corresponding eigen- 

states not presented here. The same applies to the bound 
states except for the zero energy state for ky = 0, which 
is invariant under the operation of az, and thus directly 
yields the linearly dispersing chiral state plotted as black 
straight line in Fig. [2]d. 

In the asymptotic expansion of \ipk^,ky) away from 
the domain wall, we find that finite ky only modifies the 
spinor structure and therefore the transmission coeffi- 
cients remain independent of ky. Thus, the transmission 
probability can be directly obtained from ^ by setting 

— k^ — kx, viz.. 



OOP 



(fcx) = 



sinh^(7rfcxw) 



sin (ttA) -I- si-nh {-Kk^w) 



(8) 



We observe that Tqop features oscillations in w with pe- 
riod and, in particular, for DWs with A e N, i.e. the 
reflection is completely suppressed for any kx and ky. 
The ballistic domain wall resistance reads JGqop — 
dfc,. 



J-fep 



-Roop( 



fcp 



kT) with the Fermi wave- 



2k 

vector hvkp = ^^1^— . Since the spectrum is identi- 
cal on both sides of the wall, (5Gm = here. This result 
is plotted in Fig. [3| and we recognize the oscillations in 
A originating from Tqop- For the special points where 
the domain wall width is an integer multiple of the mag- 
netic length, we find the domain wall to be completely 
transparent for the Dirac Fermions and 5G drops to zero. 

When the wall is much shorter than any other relevant 
transport length scales of the system, i.e. A, kpw <^ 1, we 
can expand the sin in Eq. ^ and perform the integration 
analytically, 



(5G, 



OOP 



2hvkpEp 



log 



Ep^ 
Ep 



hvkp 
hvkp 



(9) 



For Ep ^ M, we find the asymptotic behavior i5Goop 
{M / Ep)"^ log{2Ep / M) , while in the opposite case when 
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the Fermi level comes close to the band edge Ep « M, 
we have SGqop — >■ 1, as for the IP-configuratfon. 

In the limit of wide walls A, kpw S> 1, we obtain the 
approximate result 



SG. 



AS>1 



OOP 



log 2 

{nkpw) 



(10) 



with the modulation function J^(A) = ^J^^ which 
inherits the periodicity from Toop, in particular, J- {A) 
vanishes for integer A. We see that the envelope de- 
creases with l/w'^, similar to the result for the IP DW, 
Eq. ([t]). We find that the approximate formula fits well 
already for A > 0.5 and show in Fig. [3] the comparison 
between the full integration and the approximation (10 1. 



In a realistic system, there is always impurity scatter- 
ing which mixes channels with different k and therefore 
reduces SGm- However, as long as the wall width is much 
smaller than the mean free path, i.e. w <C /mfpj our bal- 
listic treatment is approximately correct. Impurities in- 
duce a finite scattering between the channels, however, 
due to the chiral nature of the electron dispersion, back- 
scattering is reduced. Thus, in the OOP case and when 
the DW conductance 6G is vanishing for reflectionless 
DW potentials, we expect SG to remain significantly re- 
duced in comparison to reflecting DWs for A ^ N. 

Super symmetry. The special characteristics of the do- 
main wall resistance in the OOP geometry - the periodic- 
ity in A and the perfect transmission - can be understood 
as a supersymmetry encoded in the Hamiltonian Hqop- 
Using the language of supersymmetry US] , we intro- 
duce the generalized creation — —d^ + A tanh(^) and 
annihilation operators qa — + A tanh(^) with dimen- 
sionless ^ — x/w. Introducing spin raising and lowering 
matrices a± = (cr^ ± iay)/2, the supercharges [24] are 
Qa = aAC+ and = a^(T_, so that the OOP Hamil- 



n = 

= n^-kl is di- 



tonian reads H ~ UHoopU^ — —Qa — 

We observe that the operator defined by 

{Qa,Qa} = H^Ajai} + 5[aA,a^]o"z 
agonal and essentially the square of our original Hamil- 
tonian. Solving the eigenvalue equation for H is equiva- 
lent to the one for H. Note that in this representation, 
^[aA,a|\]crz removes the zero-point energy, thus allow- 
ing for the zero energy state, see Fig. [2] Obviously, 
[H, Qa] — 0, which expresses the supersymmetry be- 
tween the two components tp^^^ of the spinor. Explicitly, 



= qaOa and Til = a^OA are supersymmetric partner 
Hamiltonians, which are iso-spectral, except that 'H^ has 
one additional bound state. Furthermore, the reflection 
and transmission coefficients defined by H^,^ differ only 
by a phase [2i] . 

For the tanh-DW profile, a second symmetry exists, 
[n - A2,CT±eT3A] = 0. This means that H^,^ are part 
of a hierarchy of form-invariant supersymmetric part- 
ner Hamiltonians, each differing from its neighbor by 
A — >■ A -I- 1. All Hamiltonians in this hierarchy have the 



same transmission-/reflection probabilities which readily 
explains the oscillations in SGqop as w/Im varies. Due 
to the scaling in ^ = x/w, there is an additional smooth 
dependence on w which yields the factor 1 /w^ in SGqqp , 
Eq. (10). Furthermore, if A G N, the constant potential 



is part of the hierarchy and thus all Hamiltonians in the 
hierarchy are reflectionless as well. Note, however, that 
the constant potential is not realized in our system, since 
the scaling becomes singular for w — > 0. Finally, the 
number of bound states differs by 1 between two neigh- 
bors in the hierarchy, which explains that the number of 
bound states is given by [AJ . We point out that similar 
reasoning can be used for the IP wall configuration, there 
however, the supersymmetric hierarchy is constructed in 
(A, A;j,)-space and thus, the fcy-integration performed in 
6G averages out the characteristic signature of the hier- 
archy. 

We remark that H^^i describes essentially a free par- 
ticle in the Poschl- Teller potential, which is known to be 
reflectionless for certain parameters. In the reflectionless 
case however, a transmitted wave still acquires a phase 
which has consequences on a wave packet passing such a 
potential: it narrows and is ahead in time as compared 
to a freely moving wave-packet [5S]. For optical systems, 
these potentials have been realized recently using arrays 
of evanescently coupled waveguides [?7] . 

Conclusions. We have analytically calculated the bal- 
listic DW conductance for in-plane and out-of-plane mag- 
netic domain walls induced into the surface states of a 
topological insulator. For the in-plane DW, the DW con- 
ductance is dominated by the spectrum mismatch im- 
posed by the opposite magnetization directions within 
the domains. For the out-of-plane DW, we unexpectedly 
flnd oscillations in the wall-width dependence with pe- 
riod Im (magnetic length). In particular, integer w/Im 
constitute a family of reflectionless potentials. We can 
understand these features using the idea of supersym- 
metry and find them to be a result of the dispersion of 
the topological surface states together with the specific 
tanh DW-profile. Detecting the oscillatory DW resis- 
tance could be a unique signature of the chiral Dirac 
surface states. This could be realized by placing an addi- 
tional top-gate on the ferromagnetic insulator which can 
be used to tune the exchange field M. The oscillations 
would still be visible if the ratio Ep/M changes slowly 
enough. 
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